Scaling in the time-dependent failure of a fiber bundle with local load sharing 
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We study the scaling behaviors of a time-dependent fiber-bundle model with local load sharing. 
Upon approaching the complete failure of the bundle, the breaking rate of fibers diverges according 
to r(t) cx (Tf — t)~ , where Tf is the lifetime of the bundle, and £ fa 1.0 is a quite universal scaling 
exponent. The average lifetime of the bundle < T/ > scales with the system size as N~ s , where S 
depends on the distribution of individual fiber as well as the breakdown rule. 

PACS numbers: 64.60.Fr, 62.20.Mk, 64.60. Ak, 05.45. +b 
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I. INTRODUCTION 



The failure of disordered materials under load is a com- 
plicated phenomenon, the modelling of which is a subject 
of great interest because it forms the basis of numerous 
applications from space technology to paper making jlj . 
The failure process also represents an important class of 
pattern formation and scaling problems The fiber- 
bundle model, as a simple and interesting theoretical 
model in this field, has been studied extensively. The 
early studies on the static fiber-bundle model might be 
traced back to the work by Daniels ||, while the time- 
dependent method to the model was proposed by Cole- 
man [0]. In a recent paper Gomez et al developed 
a probabilistic method for solving the time-dependent 
model. In the static model, each fiber in the bundle 
is assumed to have a strength threshold, a load above 
which will break it instantly, while a load below which 
does no harm. In the time-dependent model, each fiber 
is assumed to have a lifetime under a given load history, 
and it breaks because of fatigue. The load-sharing rules, 
which describe how the load of a broken element is trans- 
ferred to survival elements, are essential to the definition 
of the model. In what is called Equal Load Sharing (ELS) 
model, the total load of the bundle is equally shared by 
all survival fibers, while in the Local Load Sharing (LLS) 
model the load of a broken fiber is transferred to its near- 
est neighbors. Hierarchically organized fiber bundle was 
also proposed, and has received much attention especially 
in the geophysical literature ||[?J . Various aspects of the 
fiber-bundle model have been investigated, such as the 
strength distribution for static model and the life- 
time distribution for dynamic one QJ^]. In this paper, we 
will study an LLS time-dependent model, and investigate 
the scaling behaviors in its breaking process. 

Let us consider a fiber bundle consisting of N fibers. 
We assume that when a fiber is subjected to a load his- 
tory <j(t), some damage will accumulate, which is de- 
scribed by 



d(t) = f v[cr(T)}d: 
Jq 



(1) 



where the load-dependent v(a) is introduced as a hazard 
rate, which is usually referred to as breakdown rule Q in 
the literature. 

A fiber, say fiber i, is assumed to have an endurance 
threshold (or say, critical damage) d\, which is drawn 
from a cumulative distribution 



P(d?<d) = l-exp[-*(d)], 



(2) 



where ^(x) is the shape function. Previous theoretical 
and experimental work favors a shape function of 
the form 



*(x) = x 



- J* 



(3) 



As for the breakdown rule v(o~), two special forms are 
widely used in the literature: the power-law form 



(4) 



v p {a) = v \ — 



and the exponential form 



t> exp 



o-q 



(5) 



with vq, (To, p, </>o, t] all positive constants. 

Under load each fiber will break when the damage ac- 
cumulated exceeds its endurance threshold, and all fibers 
will break eventually, leading to the complete failure of 
the bundle. Let us denote the total load on the bundle 
by Na. In general, <r is a function of time. For exam- 
ple, it can be a linearly increasing function or a periodic 
function of time Q. In this paper, we will consider the 
simple case that a is a constant. In the following nu- 
merical calculations, if not otherwise specified, the load 
is set to be a = ctq. It should be noted that although 
the total load on the bundle is constant, the loads on the 
individual fibers <7j(£)'s are not. 



1 



II. THE LLS MODEL 

We consider a fiber-bundle model with the LLS rule. 
N fibers are arranged evenly on a circle, and each of them 
has two adjacent neighbors. The total load on the bun- 
dle No, kept constant in this study, is shared by survival 
fibers. A survival fiber i carries the load <7j = KiO, where 
the concentration factor Ki = 1 + (h + ?"i)/2. Here U ( 
r, ) is the number of broken fibers on the left (right) of 
fiber i. It is clear that J^i %i = N, so the total load 
is conserved. With such a load sharing rule, the load 
of a broken fiber is transferred to the survival neighbors 
on both sides. Note that this rule is different from the 
one-side case Jlo"[ , in which the load of a broken fiber is 
transferred only to its neighbor on one side. 

This LLS fiber-bundle model was in early years de- 
veloped by Harlow and Phoenix || to model the fail- 
ure of a unidirectional composite material under tensile 
loads. The model has ever since drawn much attention 
of many authors. In recent years, the static LLS fiber- 
bundle model was studied in terms of the burst-size dis- 
tribution |0]-[n| and the failure probability of the bundle 
under a given load 0,|l5| . In this study, we will focus on 
the scaling behaviors of this dynamic LLS fiber-bundle 
model. 



III. SCALING OF BREAKING RATE WITH 
TIME TO FAILURE 

Let Nf (t) be the number of broken fibers in the bundle 
at time t, with N f {0) = and N f (T f ) = N, where T f 
is the lifetime of the whole bundle. The breaking rate of 
the bundle is defined as 



r(f) 



St 



(6) 



We have performed extensive Monte Carlo simula- 
tions of the breaking process of the time-dependent fiber- 
bundle model with LLS, and found that in a wide range of 
parameter value the breaking rate r(t), upon approach- 
ing the complete failure, scales with the time to failure 
as 



r(t) (x (T f - ty 



(7) 



and the scaling exponent £ « 1.0 is of a quite universal 
value. Examples of the behavior of the breaking rate are 
shown in Fig. [|. In this log-log plot, dashed lines with 
slope —1 are also shown for reference. The numerical re- 
sults are not very smooth because of fluctuation, but the 
general trend of the breaking rate r(t) agrees well with 
Eq.@. 

In what follows, we try to understand the scaling be- 
havior (Q) through analytical treatment. In the discus- 
sion, we take the limit N — ► oo. Let us call the connec- 
tive broken fibers bounded by unbroken ones as a crack. 



The size of a crack is the number of broken fibers. Be- 
cause of the local load-sharing rule, the fibers bounding 
a larger crack experience heavier load than those bound- 
ing smaller ones. So when a major crack is formed in 
the bundle, breaking will mostly occur along it. In other 
words, it is the fibers adjacent to the major crack that 
will most probably break in the next step. This can be 
seen from the evolution of the size c m of the biggest crack. 
Fig. |2| shows c m versus the total number of broken fibers 
in the bundle. At the early stage of the failure process, 
c m remains constant for some time ( A ) , which indicates 
that small cracks nucleate at different locations. As more 
and more fibers break, some small cracks will coalesce or 
grow to form a major crack, and then the major crack 
grows, which is reflected in this figure by a linear increase 
of c m with Nf with slope 1(B). During its growth, the 
major crack may also coalesce with some small cracks 
and become even larger, indicated in the figure by local 
slopes steeper than 1 at some points ( e.g., C ). 

Suppose the size of the major crack is Nf(t) — k, where 
k is the number of failed fibers which do not belong to 
the major crack. The loads on the fibers adjacent to the 
major crack are [1 + (Nf — k)/2]o, so damage will accu- 
mulate in these fibers with the rate v([l + (Nf — k)/2]a). 
The breaking rate of these fibers can be assumed to be 
proportional to v(-), and one has 



it) 



dN f (t) 
dt 



A(t)v 



1 



Nf-k 



(8) 



where A(t) is a factor that depends on the accumulated 
damages in the fibers and their endurance thresholds. 
An exact calculation of A(t) is extremely difficult and 
might be impossible. We assume that the variance of 
A(t) is unimportant and take A as a constant for sim- 
plicity. The validity of this assumption is verified by the 
agreement with numerical results. Note that sometimes 
a fiber adjacent to the major crack happen to be also ad- 
jacent to a small crack, resulting in a little more load on 
it, the influence of which on the breaking rate however, 
is negligible upon approaching the complete failure. 

For the exponential form of breakdown rule (^), we 
have 



dN f (t) 
dt 



= A(j)Q exp 



V 1 + 



Nf-k 



OQ 



and therefore, 



it)=a- 1 (T f -ty 



(9) 



(10) 



where a = 7]o/(2oq), Tf is the value of time that gives 
Nf(T f ) - oo. 

For the power-law form of breakdown rules (0), 



dN f (t) 
dt 



= Au 



Nf-k 



(11) 



and 



2 



r(t) = C 



P-1 



cx(T/-t)- 



(12) 



P(T 7 < t) = 1 



-JV[i/((r)t]' 3 



(15) 



with C = Av a {a/a ) p , and N f (T f ) -> oo. So £ = 
1 + l/(p — 1). Since p is of quite large value, typically 
between 10 and 80 Q, it is not surprising that £ « 1.0 in 
the numerical simulations. 



IV. LIFETIME OF THE BUNDLE 

In deducing the scaling of the breaking rate, we have 
taken the thermodynamic limit by setting Nf(Tt) = oo. 
In numerical simulations however, we cannot realize in- 
finite system size. Given the local load-sharing rule, the 
lifetime Tf of a fiber bundle depends on the endurance of 
each fiber. Due to fluctuation, Tf is different from bundle 
to bundle. Since the fluctuation is related to the system 
size, the average lifetime < Tf > of the bundle should in 
principle depend on N, which is known as size effect. We 
found that in general the average life time < Tt > scales 
with the system size as 



< T f >oc N~ 



(13) 



where < • ■ ■ > means the ensemble average. Some of the 
numerical results are shown in Fig. ^, in which the power- 
law fit to the data is quite good. Some other forms of fit 
to the data were also tried, but none of them is better 
than the power law. It should be noted that in the static 
LLS model the average strength of the bundle follows a 
logarithmic dependence on the system size 

HOD- 

The exponent 5 for the power law, however, is not of 
a universal value. It depends on the breakdown rule as 
well as the distribution of damage endurance for individ- 
ual fiber. We performed extensive numerical simulations 
to explore the relation between the exponent 5 and the 
parameters (3, p and rj. Some results are listed in Table 
I There seems no simple general expression relating 6 to 
/3, p, and rj. For some limiting cases, however, we can get 
a simple relation. From Table |, one can see that when 
p or rj is large, the value of the exponent 5 is very close 
to 1/(3. This result can be understood by the lifetime 
distribution of the fiber bundle. When p or rj is large, 
the fiber bundle breaks in the following way: when the 
weakest fiber breaks, it will form the crack that leads to 
the failure of the whole bundle. So the lifetime of the 
bundle will be expended mostly in the weakest fiber, and 
is thus determined by it. From Eqs. (|l|), (||) and (|||), the 
lifetime of an individual fiber under a constant load a, is 
distributed as 



P(t f < t) = 1 



(14) 



For a bundle of N fibers, if the bundle's lifetime is 
determined by the lifetime of its weakest element, the 
lifetime distribution for such a bundle is, by the weakest- 
link rule and when N is large, 



And this is the Webull distribution, with which the av- 
erage lifetime of the bundle is 



< T f >- 



f-OC 

/ tdP(Tf 
Jo 



<t)= td(l-e 



(16) 



Changing the variable of integration Nt@ — t@ , one gets 



< T f >= N 



- /v-V/3 



■d(l-e- [v ^ T]f> ). (17) 



The integration in the above equation is independent of 
N, so < T f >oc N^ 1 /! 3 , and 5 = 1/0. 

From the numerical results, we notice that 5 = 1/(3 is 
not satisfied by all values of p and rj. The deviation of 5 
from 1/(3 may indicate the deviation of the lifetime dis- 
tribution from the Webull distribution. In Fig. ||, we plot 
the lifetime distribution of the fiber bundle with Webull 
axes, that is, to plot ln{— ln[l — P(t)]} versus Int. If the 
distribution is of Webull form, P(t) = 1 — exp(— at" 1 ), one 
should see a straight line in such a plot, and the slope 
of the line gives the Webull modulus m. For the case 
(3 = 2 and p = 40 (Fig. |[b), we get a quite straight line, 
and the best linear fit to the distribution curve gives the 
Webull modulus m w 2.03, very close to (3 = 2. Notice 
that for this case 5 w 0.50 = 1/(3. For the case (3 = 1 
and p = 10 (Fig. [|.a), however, the distribution curve is 
not a straight line, indicating that the lifetime is not very 
well Webull distributed. For this case 5 « 0.49, which is 
quite different from 1/(3 = 1.0. 

In the early studies on the lifetime distribution, 
Phoenix and his collaborators H were able to obtain 
an approximation to the lifetime distribution of the fiber 
bundle, which was also of Webull form. Their results 
were based on the idea that whenever a crack of criti- 
cal size, called fe*-crack in their paper, emerges in the 
system, the bundle will fail instantly. 



V. CONCLUSIONS 

In conclusion, we have studied some scaling behaviors 
of the time-dependent fiber-bundle model with LLS rule. 
In a quite wide range of parameter value, the breaking 
rate scales with the time to failure as (Tf — t) -1 - The 
average lifetime of the bundle scales with system size as 
N~ s , with 5 dependent on the breakdown rule and en- 
durance distribution of individual fiber. In the limiting 
cases that p or rj is very large, the lifetime distribution of 
the bundle can be well approximated by a Webull form, 
and the Webull modulus for this distribution is just the 
shape-function parameter (3, and the scaling exponent 
5 = 1/(3. 
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FIG. 1. The breaking rate r(t), denned in the text, scales 
with the time to failure as (T/ — , where £ ~ 1.0 is a quite 
universal value, (a) Using the power-law breakdown rule (Q) 
with p — 10, and the shape function (H) with (3 = 2. (b) 
Using the exponential breakdown rule (H) with r\ = 1.0, and 
shape-function parameter (3 = 4. In both (a) and (b), the sys- 
tem sizes are N = 100, and the dashed lines show the curves 
for y oc x^ 1 for reference. 
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FIG. 2. An example of the evolution of the biggest crack 
in the failure process of the fiber bundle. The exponential 
breakdown rule is used with rj — 1. The other parameters are 
N = 100, (3 = 4. (b) is a part of (a) enlarged. 
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FIG. 3. The average lifetime of the bundle scales with 
system size N according to a power law. The circles are re- 
sults from numerical simulations with at least 100 samples, 
the solid line is for the power-law fit y = ax~ s to the numer- 
ical data, (a) (3 = 1, n = 10, a = 8.14 x 10~ 6 , S = 0.60. (b) 
(3 = 2, p = 40, a = 0.87, S = 0.50. 
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FIG. 4. The lifetime distribution of the LLS fiber bundle. 
The results in this figure are from simulations of 10 4 samples. 

(a) j3 = 1, p = 10, N = 1000. The curve is not a straight line. 

(b) P = 2, p = 40, N = 800. The curve is a quite straight 
line, indicating a Webull distribution P(t) = 1 — exp(— at m ). 
The best linear fit to the numerical data in (b) gives the slope 
m w 2.03. 



TABLE I. The exponent 8, defined in Eq. (|T3|) , depends on 
the breakdown rule as well as the endurance distribution of 
the fibers. 





p= 10 


p = 20 


p = 40 


7,= 1 


r? = 10 


7? = 20 




1 


0.49 


0.83 


0.97 


0.17 


0.60 


0.97 


p = 


2 


0.33 


0.47 


0.50 


0.10 


0.42 


0.52 


= 


4 


0.23 


0.26 


0.27 


0.053 


0.23 


0.26 
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